Statistical Physics of the Yielding Transition in Amorphous Solids 
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The art of making structural, polymeric and metallic glasses is rapidly developing with many 
applications. A limitation to their use is their mechanical stability: under increasing external 
strain all amorphous solids respond elastically to small strains but have a finite yield stress which 
cannot be exceeded without effecting a plastic response which typically leads to mechanical failure. 
Understanding this is crucial for assessing the risk of failure of glassy materials under mechanical 
loads. Here we show that the statistics of the energy barriers AE that need to be surmounted 
changes from a probability distribution function (pdf) that goes smoothly to zero to a pdf which is 
finite at AE = 0. This fundamental change implies a dramatic transition in the mechanical stability 
properties with respect to external strain. We derive exact results for the scaling exponents that 
characterize the magnitudes of average energy and stress drops in plastic events as a function of 
system size. 



In this Letter we focus on the statistical physics of the 
yielding transition at very low temperatures and quasi- 
static external straining conditions, (the so-called athcr- 
mal quasi-static or AQS limit) where very precise simula- 
tion results exist for the dependence of energy and stress 
drops in plastic events as a function of system size [lj. 
Consider Fig. [T] which demonstrates the nature of the 
yielding transition. We plot here the conditional mean 
energy drop in a plastic event as a function of the ex- 
ternal strain 7 for two-dimensional systems (see below) 
consisting of TV particles, with TV ranging between 484 
and 20164. To read this figure properly, one should un- 
derstand that in some realizations there are no plastic 
events at all at a given external stain. What is mea- 
sured here is the size of the mean energy drop if such 
a drop happened at an external strain value between 7 
and 7 + c?7, averaged over numerous realizations of the 
random structure of the system (see below for details). 
We see that in the early stages of the loading, the plastic 
events are localized and the amount of energy released 
in events is system-size independent. This is followed by 
a smooth rise in these curves, showing an increasingly 
sharper transition to the plastic flow state in which the 
plastic events become non-localized avalanches whose to- 
tal energy release increases with the system size. This 
very interesting system size dependence will be quantified 
below. We note in passing that the stress itself cannot 
be a proper order parameter; states with the same stress 
level (shown for example in Fig. [T]as two magenta circles) 
have very different conditional mean plastic energy drops. 
Here we explore the statistical physics that is responsible 
for the difference between these iso-stress states, which 
also have very similar potential energy and pressure. We 
point out that the precise nature of this strain-induced 
transition from the solid-like jammed state to the steady 
flow state, where the plastic flow events resemble liquid- 
like dynamics, is still unclear. Although the increasing 
availability of computational power has recently led to 
many important observations and conclusions regarding 
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FIG. 1: Color online: Evolution of the conditional mean en- 
ergy drops with the loading, for all system sizes simulated, 
increasing from bottom to top. Superimposed (scale on the 
right ordinate) is the mean stress vs. strain curve for the 
largest system of N = 20164. The magenta dots represent 
equi-states states but of highly different stability properties. 



the statistics of the steady flow state [IrEl], a clear-cut 
identification of the physics that control the approach to- 
wards steady state, has not been presented yet. The aim 
of this Letter is to close this gap and to offer some exact 
results. We stress that this desired analysis is best con- 
ducted in the AQS limit since much is known there about 
the nature of the plastic events themselves, as these are 
determined by mechanical instabilities which can be seen 
as a saddle node bifurcation in which the lowest eigen- 
value of the Hessian matrix going through zero 0B| • De- 
note the potential energy of the system as U(ri) where r, 
are the positions of the particles, and the Hessian matrix 
as Hij = d 2 U / dridrj. The Hessian is a real symmetric 
matrix; we denote its lowest eigenvalue (excluding the 
Goldstone modes) as Xp . It was established |4|, |6| that 
when the external strain 7 reaches a critical value jp, 
Xp vanishes with a square root singularity, i.e. such that 
Xp oc y/^yp — 7. We will show that this simple singularity 
determines the numerical values of a number of interest- 
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ing exponents that appear in the statistical analysis. 

Below we employ a model glass-forming system with 
point particles of two 'sizes' but of equal mass m in two 
and three dimensions (2D and 3D respectively), inter- 
acting via a pairwise potential which had been fully de- 
scribed in [8(. The experiments performed are as fol- 
lows: For un-deformcd isotropic systems we measured the 
strain at which the first plastic event takes place, and de- 
noted it as A7; so . Each such measurement was performed 
on a freshly produced amorphous solid, quenched from 
the high temperature liquid at the rate of 5 x 10 -5 ^:. 
Then the AQS scheme (see Q for details) was utilized to 
strain the system up to the first mechanical instability 
occurring at some strain value A7i so . Statistics of A7i so 
were collected for a variety of system sizes, see below. 
In the elasto-plastic steady state we first strained sta- 
tistically independent systems for 100% strain to reach 
stationarity, and then collected statistics as shown below. 

In the steady flow state, the statistics of the energy 
drops AU, the stress drops Ac and the strain intervals 
between successive flow events A7 become stationary. 
Quite surprisingly, one finds that the averages of these 
quantities obey scaling relations with the same exponents 
in two and three dimensions: 

(AU) ~ lN a , (Act) ~ sN? (1) 
(A 7 ) ~ . (2) 

In Fig.[5]the mean energy drop (AU) and mean strain in- 
terval (A7) for our model system are displayed, together 
with the scaling laws ([1]). In the upper panels we show 
results in two dimensions and in the lower panel in three 
dimensions, and in both a s» 1/3 and /3 ps —2/3. A scal- 
ing relations a — /3 = 1 was already established before 
Q. In this Letter we propose that the respective values 
1/3 and -2/3 are exact. 

The yielding transition is underlined by the fact that 
for the first plastic event when strained from a freshly 
quenched isotropic state the statistics is entirely different, 
with AU ~ ./V representing a localized event without 
any system size dependence. On the other hand the first 
plastic event does not occur for any infinitesimal value of 
7 and careful measurements of the mean strain interval 
(A7 iso ) that separates the un-deformed state from the 
first plastic event results in a scaling law 

A 7iso ~ N&*° , /3 iso « -0.62 . (3) 

How can we understand the difference between /? and /?i so 
and what determines their numerical values? 

Starting from any given mechanically stable state, at 
the un-deformed state or at the steady state, the system 
has a set of O(N) energy barriers AE which are coupled 
to the external strain. One of those needs to be sur- 
mounted in order to have a plastic event. In AQS condi- 
tions the one chosen will be the one which has the small- 
est A7i so (in the isotropic state) or the smallest 7p — 7 
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FIG. 2: Color online: panel (a). Mean energy drop (At/) 
and mean strain interval (A7} in 2-dimensions as functions 
of system size, measured in AQS simulations of steady plastic 
flow of a model glass former, see text. Panel(b): the same for 
three dimensions. The continuous lines represent the scaling 
laws {!]) and ([2]). The scaling exponents are the same in 2D 
and 3D. 
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FIG. 3: Color online: Scaling of energy barriers for various 
system sizes, see for details. The slope of the continuous 
lines is 3/2. 

(in the steady state) . As a function of the external strain 
this barrier reduces until it vanishes at the saddle node 
bifurcation where Ap vanishes In Refs. 0, 0] it 

was shown that the manner in which the energy barrier 
vanishes is determined by the saddle node singularity. In 
other words, it was established that close to 7p 

AE cx A3 ~ ( 7 p - 7 ) 3 / 2 . (4) 

We stress that this result is valid, sufficiently close to 
7p, equally well when starting from equilibrium, where 
7p represent the value of the strain for the first plastic 
event, or in the steady state, where 7p is any value of the 
strain where a plastic event occurs. It turns out that the 
scaling law ((4]) is obeyed, at least in the class of models 
in which the potential is purely repulsive, for a very long 
range of 7p — 7, see Fig. [3] We will use this in this Letter, 
coming back to the question of universality at the end. 
In terms of distributions, the possible plastic events 
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can occur anywhere in the system, and the number of 
possible sites increases linearly with A. However, we re- 
alize that every time that we observe a plastic event in 
AQS conditions, it is one of the lowest barriers out of 
0(A) barriers that is overcome [10(. We cannot directly 
measure the distribution of energy barriers but need to 
concentrate on the extreme statistics regarding the min- 
imal value. In the isotopic state, we expect the distribu- 
tion of putative strain values x that could be assigned to 
mechanical instabilities to be Poissonian, of the form 



p(x) = x v h{x)/Z 



Z : 



x v h(x) T) > . (5) 



where h(x) decays rapidly for i» 1. It is well known 
then that if we now take a set of A 3> 1 independent 
samples from such a distribution, then the probability 
distribution g(y,N) of the minimal element of the set 
(denoted y = A7 iso ), is the Weibull distribution [11 1 
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In this equation yo ~ ]y- 1 /( 1 +">i) j s ^he mean value of y 
with respect to the Weibull distribution. A test of this 
logic is presented in Fig. [4] for systems in 3-dimensions 
where we show the distribution of A7j so of isotropic 
states, in excellent agreement with Eq. ([5]) with r\ w 0.6. 
In 2-dimensions the scaling with r\ 0.6 is also recov- 
ered. We expect the distribution of the coefficients of 
the scaling law Eq. ((4]) to be regular, implying that the 
probability distribution of the energy barriers has, for 
low values of AE, the form 



p(AE) ~ (AE)* , 77 = (27 7 - l)/3 



(7) 



For rj w 0.6 we find f] ~ 1/15 > 0. This is consis- 
tent with the notion of solidity; one expects that for a 
solid the probability of finding a zero barrier is strictly 
zero. Finally, the mean value of the minimal values of 
AE scales as predicted by the Weibull distribution, i.e. 
(AE) ~ JV-VCl-H). Substituting this in the mean of 
Eq. dU) we recover Eq. ([3]) with the observed exponent 
fii SO s» 0.62. We thus conclude that the first plastic event 
is dominated by extreme statistics of the minimal barri- 
ers for plasticity with a pdf of the energy barriers that 
goes to zero for AE — > as is expected from a solid. 

The picture changes dramatically when we investi- 
gate the elasto-plastic steady state. There the distri- 
bution of strain intervals between successive events A7 
changes qualitatively as seen in Fig. H Now the power- 
law tail of the distribution is too shallow to imply that 
limAB->o p( AE) = 0. Indeed, we propose that the hall- 
mark of the yielding transition is that in the elasto-plastic 
steady state the probability to find a zero value of the 
energy barrier is non-zero, cf. 12j . In other words, the 



criterion of solidity is no longer applicable. We also can- 
not expect that the statistics of A7 follows the Weibull 




FIG. 4: Color online: the measured pdf (in 2-dimension for 
N — 4096 in the panel a) and in 3-dimensions for N = 2000 
in panel b) for the observed value of x = A7i so for isotropic 
systems (in red), and of x = A7 for the steady state (in blue). 
The scaling law for the mean shown in the inset indicates 
a value of 77 ~ 0.6 for both 2 and 3 dimensions. Note the 
dramatic change in the power-law tail of the distributions at 
small values of A7: in equilibrium the tail guarantees that 
the probability to see a zero value of AE is zero. This is not 
the case in the steady state, and this is the physical hallmark 
of the yielding transition. The black line through the red dots 
in panel b) is the Weibull distribution (J(j]). 



distribution, since the values of 7p in subsequent plastic 
events become highly correlated and history dependent. 
Indeed the steady state distribution shown in Fig. [4] can- 
not be fitted to a Weibull distribution. On the other hand 
we can still expect that the numerical values of the min- 
imal energy barriers that need to be surmounted remain 
statistically independent due to the avalanches, and see 
next paragraph. The consequence of this is that with A 
independent random sampling of AE with a finite prob- 
ability to find AE = 0, the scaling of the minimal value 
must scale like 1/N. Using this in Eq. (j4j leads to the 
proposed exact values of a = 1/3 and j3 = —2/3. 

To understand the pdf of strain intervals as measured 
in the steady state (cf. Fig. 0J, we need to explain how 
the avalanches that occur after every plastic event re- 
fresh the statistics and renormalizes the pdf to the form 
seen in Fig. [4j The cumulative energy associated with 
an avalanche grows sub-extensively with the system size, 
(like A 1 / 3 ) and therefore the impact of these avalanches 
remains pertinent in the thermodynamic limit. To model 
the effect of the avalanches imagine that we consider a 
population of A energy barriers AEi sampled from a dis- 
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FIG. 5: Color online: pdf of the observed strain intervals be- 
tween avalanches compared to the converged pdf of the itera- 
tive model, see text for details. Here q = 2/3. The agreement 
indicates the robustness of the model and the crucial role of 
the avalanches in renormalizing the pdf. Here we opted to 
show the pdf of strain intervals because it is the directly mea- 
surable quantity, related to the energy drop by Eq. Q. 

tribution p(AE) which satisfies the two conditions: (i) 
p(AE) ~ AE° for AE — > and (ii) p(AE) -> at least 
exponentially for AE 1. The series of plastic drops 
is then modeled by the following iterative steps; at each 
step we repeat the following operations: 1) Find AE m m 
and record it. 2) For every AEi =^ A£" m i n transform, 

A25i<- (AE* /3 - AE% 3 D ) 3/2 . (8) 

This step takes into account the fact that what is changed 
in the simulations is the external strain rather than the 
energy barriers, and we have used Eq. 3) Remove 

A-E m i n and re-assign a new number from p(AE) instead 
of it. 4) Randomly remove qN numbers from the set 
(q G (0,1) is some fraction), and re-assign new num- 
bers instead of them from p(AE). This iteration scheme 
is readily performed numerically, leading to a converged 
pdf that is shown in Fig. [5] in comparison to the simula- 
tional pdf. We find an excellent agreement which under- 
lines the crucial effect of avalanches in partially destroy- 
ing the correlation between subsequent values of AE mul . 
Importantly, the scaling of the mean value is invariant 
to the choice of q in the large N limit. It should be 
noted that when the same iteration procedure is per- 
formed with a pdf that goes to zero at zero like AE V (even 
for very small 77), the resulting converged pdf is qualita- 
tively different, preserving the scaling of the mean value 
(AE m i n ) ~ ]y- 1 /( 1 + r i) _ xhis stresses the importance of 
the fundamental physics of the yielding transition which 
takes the system from a solid to liquid-like state with a 
qualitative change in limA.E->.oP(A-E). 

The main point of this Letter is that the yielding tran- 
sition is characterized by a qualitative change in the na- 
ture of the pdf 's of the energy barriers for a plastic event 
in the AQS limit. For the solid the probability to see a 
vanishing energy barrier is zero. For the elasto-plastic 



steady state this probability is finite. Physically, this 
transition is the reason for the avalanches that are ob- 
served in the steady state - there are many localized con- 
figurations with close to zero energy barrier to surmount, 
and any plastic drop anywhere in the system will find it 
easy to cause all these configurations to cross the instabil- 
ity threshold. The implication of this qualitative change 
is the change in scaling exponents which are shown in 
Eqs. ([2j> and ((3]). One main result of the Letter is the 
derivation of the exact values of the exponents a and /3 
in Eq. We stress that this exact result rests entirely 
on the availability of the scaling law ((4]) . Clearly this law 
is asymptotically true for 7 — > 7p, equally well in 2D 
and 3D, implying that a and /3 are dimension indepen- 
dent. However the range of strain values over which this 
law pertains in the present study is quite amazing, and 
there is no guarantee that this will remain true in other 
models which have attractive terms in the potential of 
molecular degrees of freedom. It is therefore worthwhile 
to continue to explore the scaling properties of different 
models in the AQS conditions to delineate the existence 
of universality classes. 
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